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Abstract 

By considering a correlated two-electron transport process (TET) and using 
a diagrammatic analysis within the Keldysh nonequilibrium Green's function 
formalism, we discuss a novel charge pump by which carriers are pumped 
from a contact with low chemical potential to another contact with a higher 
potential. The TET process involves two correlated incident electrons scat- 
tering and exchanging energy with each other. The process can significantly 
affect charge current density and it involves high empty states and/or low 

filled states of the Fermi liquid of the leads. 
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In the past two decades, quantum electron transport through meso- and nanoscopoic 
structures have received considerable attention [|lj] . Charge transport theories for nanostruc- 
tures typically consider the situation of an incident electron from a lead that is scattered 
in the device scattering region and transmitted to other leads or reflected to where it came 
from. The scattering region can be rather complicated and may involve electron-electron, 
electron-phonon, electron-impurity, and other interactions and scattering mechanisms. The 
scattering region itself can be a semiconductor nanostructure such as a quantum dot (QD), 
a carbon nanotube, or a single molecule. The scattering processes may involve exchange 
of energy or spin. While situations and physics can vary in a wide range of ways, charge 
transport in nanostructures has largely, so far, been considered as involving a single elec- 
tron incident from a non-interacting lead of the device, it traverses through a complicated 
scattering region in which interactions occur, and finally exits the device through the non- 
interacting leads. We will refer them as single electron tunneling or single electron transport 
(SET) processes. 

In this paper, we go one step further by investigating correlated two-electron tunneling 
and transport (TET) processes which, to the best of our knowledge, has not been studied 
before. In this process, two incident electrons with energies e\ and e 2 from the leads tunnel 
into the scattering region of a device, they scatter with each other and exchange energy 
(they could also scatter and exchange energy with other particles in the scattering region). 
Following the scattering event, they transmit to the outside world through the contacts, 
but with energies e x and e 2 which are different from the initial values e±, e 2 due to energy 
exchange. If e\ + e 2 = e x + e 2 , i.e. the two incident electrons do not exchange energy with 
other particles, the transport process as a whole can be considered as elastic. In the TET 
process, the variation of electron number in the scattering region is always two: a situation 
that is different from SET. We show that such TET processes can indeed occur, and it 
affects measurable physical quantities such as charge current. In particular, TET induces a 
transport vertical flow — meaning carriers incident with energy e\ exit at a different energy 
e[, resulting to a non-conserved current density, i.e. 2nJn( e ) 7^ 0. Using the characteristics 
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of the correlated two-electron transport process, we design an interesting charge pump so 
that incident carriers can be pumped from a contact with low chemical potential to another 
contact having a higher potential. 

We consider a QD coupled to two or three leads, described by the following Hamiltonian: 

[tnkai ka d a + H.c.}, where a} nka (a nka ) 
and d^ a (d a ) are creation (annihilation) operators in lead n and in QD, respectively. The QD 
includes two states and has an intradot Coulomb interaction U. Subscript a is the spin 
index, it may also indicate other quantum numbers. To account for a magnetic field, we let 

In the following we focus on investigating the elastic TET by analyzing the behavior 
of current density j n (e) from lead n to the QD. j n (e) relates to the current I n through 
In — I jn(t)de. We also define the electron occupation number density operator N n (e,r) = 
J2k,a I eteta \ika{ T ) a nka{T + ^)|^- N n (e, r)de describes the electron occupation number in lead 
n in the energy range e to e + de at time r. Current density j n (e) can be calculated from the 
time evolution of N n (e, r): j n (e) = — e < J^iV n (e, r) >. This leads to (in units of h — 1): 

j n (e) = -elmY^- Wn{e)G r a {e) + G<(e)] (1) 

where linewidth function T n (e) = 2tc J2k |^nfc| 2( 5(e — e nk)', /n(e) is the Fermi distribution 
function for lead n; G r ^ < {e) are the retarded and Keldysh Green's functions of the QD. 
Using the standard equation of motion technique, G r a (e) has already been solved in previous 
work: |i 

1 + UA a n £ 
e-e Q -£i 0) (e) + £M Q (£i + ^) 

where A a (e) = [e - e a - U - E^(e) - E& } (e) - E^e)]" 1 ; E^ is the lowest-order self- 
energy, Y>^'^' a ' b are the higher-order self-energies; || and n a = Im J |^G^(e) is the intradot 
electron occupation number of state a. It is worth to mention that if temperature T is lower 
than the Kondo temperature T%, the solution of Eq.(2) has a Kondo resonance at the Fermi 
level which was the subject of many previous studies, pfl 



GKe) = - (2) 



We investigate TET processes at a temperature higher than Tk- To this end we need 
to solve the the Keldysh Green's function G^(e). Note that if one applies the commonly 
used ansatz for interacting lesser and greater self-energies 0, or using the large-t/ limit 
non-crossing approximation 0|J to solve G„(e), the two-electron scattering will be lost in 
these approximations. Therefore a more precise analysis is needed in our problem and we 
proceed as follows. Introducing the intradot electron occupation number density operator 
N a (e, t) = J e le *c?] 1 ,(r)(i Q ,(r + t)^- and due to the steady state condition < jpN a (e, r) >= 0, 
we have: 

-e£ ^ [Gt + UG r a - Gl)\ = i a (e) (3) 

where i a (t) gives the intradot verti- 

cal flow: i a {e) = ieUfe™ {< S a (0)4 (0)d a (0) d a (t) >-< S a (0)S s (t)d & (t)d a (t) >} f . The 
quantity i a (e) can be viewed as an intradot "vertical" current density jl| at energy e that is 
contributed by carriers with other energies due to the e-e scattering process. Summing over 
the index a, Eq.(3) reduces to J2 n 3n( e ) + Z)a^a(e) = 0, which is exactly the steady state 
current conservation equation so that the total current density — including the vertical flow, 
through the intradot energy level e is zero. If G r a (e) and i a (e) have been solved, from Eq.(3) 
G<(e) can be obtained immediately. 

Next, we solve the vertical flow quantity i a (e). We introduce two particle contour- 
ordered Green's functions B a {t, 0) and B c (t, 0), B a (t, 0) = - < T c [d a (t)d s (0 + )4(0)4(0)] > 
and B c (t,0) = — < Tc[d a (t)da(t)dlt(0)da(Q)] >. Although there are four operators in the 
definition B a / c (t,0), only two time indices at (t, 0) appear. We can therefore write the 
contour ordered quantities as B^f(t, 0), B^Z, B~t, and B~7. Their Fourier transformations 
can be defined as B(e) = J e iet B{t, 0)dt. Using the two-particle Green's function, the vertical 
flow i a (e) is reduced to i a (e) = ——ImB£~(e). 

The contour Green's function B a / C is solved by means of a Feynman diagram ex- 
pansion using the Wick's theorem. We take the interacting part of the Hamiltonian as 
Hi = Udtd-fdidi, and the remaining part H — Hi as the non-interacting Hamiltonian Hq. 
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The first-order irreducible self-energy which we consider is shown in Fig.l. It is clear that 
the first-order graph describes a two-particle propagation involving the exchange of an in- 
teracting energy U. Two important points should be mentioned here, (i) The vertex we 
calculate sums up all reducible diagrams constructed by the irreducible self-energy of Fig.l 
(also see Eq.(4)). This level of approximation is equivalent to that of a typical random 
phase approximation, i.e. we compute the irreducible self-energy up to order It is 

reasonable to neglect other higher order irreducible diagrams which are of orders U~ 2 and 
higher, because the interaction energy U is large, (ii) Single solid lines in Fig.l stand for the 
intradot contour-ordered Green's functions G a of the Hamiltonian H (not H ). This means 
that we have summed over all terms in the Feynman diagram expansion of G a . 
The corresponding equation for Fig.l is 



Finally, we get B+~(t,0), B~~(t,0), and B+~(t,0) from Eq.(4), and upon taking a Fourier 
transformation, £>+~(e) can be expressed in terms of G+ + (e), G^(e), G~ + (e), and G~~(e)— 
these four Green's functions are directly related to G r a and G<- This completes the analytical 
derivations. 

From Eqs. (2,3,4), the intradot occupation number n a and the vertical flow i a (e) are 
determined self-consistently. Let's consider a two-probe (n = L, R) device and the wideband 
approximation in which T n (e) is independent of energy e. Fig. 2 plots the vertical flow 
i(e) = J2 a ia(t) (solid line) as well as the current density ]l/r{.^) at a high bias. Clearly, the 
vertical flow i(e) is non-zero due to the TET processes, and its value |i(e)| has four peaks 
at energies e-f, e|, e-f + U, and e± + U, respectively. This means that the incident electron 
indeed can vary its energy in QD by e-e scattering. In contrast, in a typical SET process 
the electron keeps its energy and does not induce any vertical flow i(e). 

The physics of the TET process that induces the vertical flow is shown by inset (a) of 
Fig. 2. To start, two incident electrons from left lead having energies e± and e-f + U tunnel 



B c / a (t, 0) 



G a (t,0)G dl (t/0 + ,0) 




(4) 
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into the QD. They scatter with each other inside the QD and exchange energy to final states 
e T and ej, + U. Afterwards they tunnel out of the QD. In the vertical flow curve of Fig.2, 
two peaks at and e^ + U are negative (dips), and the other two peaks at e-f and e± + U are 
positive: precisely indicating the transfer of states from the initial ones at e± and e-j- + U to 
the final ones at e-j- and t^ + U . It is also worth mentioning that besides the new TET process, 
the usual SET processes also exist in charge transport through the QD in the present case. 

So far we have demonstrated that TET processes can exist. In the following we analyze 
several important questions concerning TET. What is its consequence? (i) TET makes 
current density a non-conserved quantity, i.e. J2 n jn( e ) 7^ 0. Of course, the total current 
is still conserved, i.e. J2 n In = 0. This can be easily proved from the definition of the 
vertical flow i(e), namely i(e) has the property J i(e)de = 0. (ii) TET can involve high 
energy empty states, namely states which are higher by about U than the highest chemical 
potential max([/,L, Hr) (see inset (a) of Fig.2); it may also involve electrons deep inside the 
Fermi sea, namely states which are lower by about U than the lowest chemical potential 
min(fj,L, /ir) (see TET process shown in inset (b) of Fig.2). (iii) TET may induce a current 
density that is flowing out from the high voltage terminal, i. e. the left lead (indicated by the 
negative peak at e± + U in jL(e) curve, and by the arrow A in inset (a) of Fig.2). Similarly, 
TET may also induce a current density that is flowing in from the low voltage terminal, i.e. 
the right lead (see arrow A in inset (b) of Fig.2). These characteristics are rather different 
from the typical elastic SET processes. 

Under what conditions does TET or the vertical flow i(e) exist? (i) We found that 
an increase (decrease) of temperature T or linewidth Y will widen (narrow) the peaks of 
vertical flow, but does not affect peak positions and heights significantly, (ii) If U = 0, 
the vertical flow i(e) = identically: TET crucially depends on this parameter, (iii) If 
U — > oo, i(e) tends to zero. This is because at large U, the intradot two-electron occupation 
is prohibited therefore TET is blockaded. In this case only SET processes occur, (iv) When 
bias potential eV = /il — is less than U, i(e) decreases drastically. In the limit of eV = 0, 
= Jn(e) = 0. 
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Are observable quantities of charge transport affected by the TET process? (i) Clearly 
the current density Jl/r^) is affected significantly as already discussed above, (ii) In general, 
the current, conductance, and n a will be affected significantly by TET (see below). The 
current noise, which reflects the e-e time correlation, will increase due to TET processes. 
However, if one uses the wideband approximation, the TET dependence in charge current 
will be lost. 

In the rest of this paper, we apply the property of TET to design a device so that 
electrons can be pumped from a lead with lower chemical potential to another lead having 
a higher chemical potential. Consider a device with three leads (n = 1, 2, 3) and consider 
the non-wideband case. We use a model of quasi-square bands where the coupling T n (e) = 
T/{exp[(\e — c n \ — W)/0.05] + I}, the width of the band is set by 2W = 1 and its center 
at c n which is dependent on the terminal voltage eV n = fi n but fi n — c n is kept fixed. More 
specifically, let's assume lead I to be a p-type semiconductor with fj,± — c\ = 0.4; lead 2 
an n-type semiconductor with /i 2 — C2 = —0.4; and lead 3 a metal with /i 3 — C3 = (see 
Fig. 4). The energy diagram of the device is set by external voltages as that shown in Fig. 4 
so that /i2 > A*i > ^3- The current density j n (e) in this case is shown in Fig.3. We note 
that ji(e) (dotted line) has two positive peaks at e± and + U; J2(e) and .73(e) each has 
one negative peak at ej_ + U and ej, respectively; and i(e) = — Z)nJn(e) has two negative 
peaks at e± and ej + U, and two positive peaks at + U and e-f. The current I n = J jn(t)de 
is quite large. We emphasize two points for this pump, (i) In this device the SET process 
almost does not occur because bands of different leads do not overlap. Then, clearly, the 
large current I n originates from the TET process: /„ — * if the vertical flow i(e) — > 0. This 
demonstrates that TET can significantly affect charge current in the general case of non- 
wideband coupling, (ii) The charge current in the terminal with the highest bias voltage, 
e.g. lead 2, is negative (solid line in the inset of Fig.3), which demonstrates the pump effect. 
The pumps works because when an electron tunnels from lead 1 to 3, it emits energy U to 
pump another electron from lead f to lead 2, through the TET process. 

More clearly, the working principle of the TET pump is summarized in Fig.4. (a) We 
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start from the situation where no charge is in the QD so that levels e-f and are empty. 
In this situation, an electron in the Fermi sea of lead 1 having energy e_f can easily tunnel 
into the QD (Fig. 4a). (b) After this electron tunnels into the QD and occupies the QD 
level of e|, the other intradot level e-f is raised up to e-f + U, so that another electron with 
energy e-f + U in lead 1 tunnels into the QD (Fig. 4b). (c) When the second electron comes 
into QD, due to the e-e Coulomb interaction U, the level e_f with its electron is raised up 
to ej_ + U, leading to a negative peak at e| and a positive peak at ej_ + U in the vertical 
flow curve i(e). Now the intradot two-electron system has total energy e-f + ej_ + U (Fig. 4c). 
Afterwards the first electron in state ej_ + U easily tunnels to lead 2 and takes away energy 
ej_ + U. The net effect is that the two electrons exchanged energy U, which is the TET 
process discussed above. When the first electron leaves the QD, the other electron at e-f + U 
falls down to e-f, leading to a negative peak at e-f + U and a positive peak at e-f in the curve of 
i(e). (d) Finally, the second electron tunnels to lead 3 (Fig.4d) and our device returns to its 
initial conformation of Fig.4a. This way an electron is pumped from lead 1 to lead 2, where 
Hi < fi2, via the TET process. We emphasize that each tunneling event from Fig.4a-d is a 
first-order normal tunneling event in which tunneling occurs at two aligning states (not like 
higher-order virtual co-tunneling process) J7|. We therefore conclude that the TET process 
should have large probability to occur so that i(e) can be near the unit value e/h (see Fig. 3). 
We have also investigated the terminal voltage (e.g. V3) dependence of current /„, shown 
in the inset of Fig.3. As V3 is increased so that eV 3 = fi 3 passes the lowest resonance state 
e-f, the tunneling event in Fig.4d can not occur and the TET process is blockaded, leading 
to a significant reduction of all currents I n (including the pumping current 1%) (see inset of 
Fig.3). 

In summary, we have investigated the two-electron correlated scattering process in meso- 
scopic system. TET induces a vertical flow in the scattering region so that electrons enter 
and exit the device with different energies. TET is found to affect current density signifi- 
cantly, and the process can involve high empty states or/and low filled states of the leads. 
The properties of TET suggests an interesting working principle of an electron pump which 
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pumps charge carries from a lead with low chemical potential to another lead with a higher 
chemical potential. In fact, if the bands of lead 1 is full and bands of leads 2 and 3 are 
empty, these results are not affected. Our proposed pump is very different as compared to 
the electron-photon or parametric pumps ||||. It should be experimentally feasible even 
for devices fabricated in two-dimensional electron gas. In that case, the bands in our theory 
can be replaced with Landau levels. 
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FIGURES 

FIG. 1. Relevant Feynman diagrams considered in this work. The single solid lines, doubles 
solid lines, and the wave lines stand for the Green's functions G a , B a / C , and interaction U, respec- 
tively. 

FIG. 2. jL/ij(e) and i(e) versus energy e in the wideband limit. The parameters are: e-f = —0.5, 
e_L = 0.1, U = 1, T = T = 0.1, and jil = — /xr = 0.7. Insets (a) and (b) are schematic plots for two 
kinds of TET processes. 

FIG. 3. The main plot shows j ra (e) versus e with ^3 = —1.3 and the inset shows I n versus 
e ^3 = fi3 for the non-wideband case. Other parameters are: e-\ = —1.1, = —0.1, U = 1.2, 
T = T = 0.05, /Lti = 0.4, and ^2 = 0.9. The dotted, solid, dashed curves in the main plot 
correspond to ji(e), ,72 (e), and j3(e), respectively; they correspond to Ii,l2,h in the inset. 

FIG. 4. Schematic plots for the working principle of the TET charge pump. 
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